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ABSTRACT 

We report on 8.4 GHz VLBI observations of SN 2001gd in the spiral galaxy NGC 5033 
made on 26 June 2002 (2002.48) and 8 April 2003 (2003.27). We used the inter- 
ferometric visibility data to estimate angular diameter sizes for the supernova by 
model fitting. Our data nominally suggests a relatively strong deceleration for the 
expansion of SN 200fgd, but we cannot dismiss the possibility of a free super- 
nova expansion. From our VLBI observations on 8 April 2003, we inferred a mini- 
mum total energy in rclativistic particles and magnetic fields in the supernova shell 
of -E m in=(0.3-14) x 10 47 ergs, and a corresponding equipartition average magnetic 
field of -E> m i rl = (50-350) mG. We also present multiwavelength VLA measurements of 
SN 2001gd made at our second VLBI epoch at the frequencies of 1.4, 4.9, 8.4, 15.0, 22.5, 
and 43.3 GHz. The VLA data are well fit by an optically thin, synchrotron spectrum 
(a = —1.0 ± 0.1; S v oc v a ), partially absorbed by thermal plasma. We obtain a su- 
pernova flux density of (1.02±0.05) mJy at the observing frequency of 8.4 GHz for the 
second epoch, which results in an isotropic radio luminosity of (6.0±0.3) x 10 36 ergs s _1 
between 1.4 and 43.3 GHz, at an adopted distance of 13.1 Mpc. Finally, we report on 
an XMM-Newton X-ray detection of SN 2001gd on 18 December 2002. The supernova 
X-ray spectrum is consistent with optically thin emission from a soft component (as- 
sociated with emission from the reverse shock) at a temperature around 1 keV. The 
observed flux corresponds to an isotropic X-ray luminosity of Lx = (1-4 ± 0.4) x 10 39 
ergs s _1 in the 0.3-5 keV band. We suggest that both radio and X-ray observations 
of SN 200 lgd indicate that a circumstellar interaction similar to that displayed by 
SN 1993J in M 81 is taking place. 

Key words: galaxies: individual: NGC 5033 - radio continuum: stars - supernovae: 
individual: SN2001gd - X-rays: individual SN 2001gd 



1 INTRODUCTION 

SN2001gd in NGC 5033 was discovered by iNakano et alJ 
(2001) on 24.82 November 2001; its explosion date is un- 
certain. The supernova had a visual magnitude then of 
14.5, and was located ^3' north-n orthwest of the nucleus 
of NGC 5033. INakano et ail fcOOll) reported the following 
position for SN 2001gd: a=13hl3m23s.89, <5=+36°38T7" .7 
(equinox J2000.0). They also reported that there was no 
star visible at the above position on earlier frames taken be- 
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tween 1996 and April 2001. An optic al spectrum obtained b y 
P. Berlind on 4.52 December 2001 dMatheson et all 1200 llh 
showed SN 2001 gd to be a Type lib s upernova well past 
maximum light. iMatheson et alJ i200ll) pointed out that 
the spectrum was almost identical t o that of SN 1993J ob - 
tained on day 93 after explosion jMatheson et alJl200d) . 
Since SN 1993J was a strong radio emitter, it was natural 
to expect SN 2001gd w ould also be a strong radio emitter. 
IStockdale et alJ J2002t) detected SN 2001gd on 8 February 
2002 at cm- wavelengths with the Very Large Array (VLA). 
Their continual monitoring of SN 2001gd since its first radio 
detection has confirmed that SN 2001gd is similar in its ra- 
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dio properties to SN 1993J (peak luminosity at 6 cm, time 
to 6 cm peak, spectral index; Stockdale et al l2003h . as in its 
optical properties. Given the similarities between the optical 
spectrum of SN 2001gd on 4.52 December 2001 and th at of 
SN 1993J at an age of 93 days jMatheson et al.ll200ll) . we 
assume throughout our paper that the supernova exploded 
93 days before 5 December, therefore fixing that event at 
t = 3 September 2001. 

A preliminary report on our radio findings is presented 
in IPerez-Torres et al] ]2005h . Here, we present the results 
of our two-epoch VLBI observing campaign on SN 2001gd, 
complemented with VLA observations and XMM-Newton 
archival X-ray data. The paper is organized as follows: we 
report on our VLBI and VLA radio observations, and on 
XMM-Newton archival X-ray data in Sect. [^1 we present our 
VLBI images of SN 2001gd in Sect. [3] we present and discuss 
our results in Sect. 2] We summarize our main conclusions 
in Sect. |S] We assume throughout our paper a distance of 
13.1 Mpc to the host galaxy of SN 2001gd, NG C 5033, based 
on its redshift (z = 0.002839: lFalco et al.ll99sl^ and assumed 
values of Ho = 65 km s _1 Mpc -1 and qo=0. At this distance, 
1 mas corresponds to a linear size of 0.063 pc. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 VLBI measurements 

We observed SN 2001gd on 26 June 2002 (2002.48) at a 
frequency of 8.4 GHz, using a VLBI array that included 
the following 12 antennas (diameter, location): The VLB A 
(25 m, 10 identical antennas across the USA), Green Bank 
(100 m, WV, USA), and Effelsberg (100 m, Germany). We 
also observed SN 2001gd on 8 April 2003 (2003.27) at the 
same frequency with the above array, complemented with 
the phased- VLA (130 m-equivalent, NM, USA), and Medic- 
ina and Noto (32 m each, Italy). However, for our second 
epoch we had to discard three antennas (Hancock, North 
Liberty, and Noto) due to the poor quality of their data. 
The duration of each observing run was eight hours. The 
telescope systems recorded right-hand circular polarization 
(RCP), and used a bandwidth of 64 MHz, except for the 
VLA (bandwidth of 50 MHz). The data were correlated at 
the VLBA Correlator of the National Radio Astronomy Ob- 
servatory (NRAO) in Socorro (NM, USA). The correlator 
used an averaging time of 2 s. 

Since previous VLA observations of SN 2001gd showed 
it to be fainter than 6 mjy at 8.4 GHz, our VLBI observa- 
tions were carried out in phase-reference mode. SN 2001gd 
and the nearby International Celesti al Reference Frame 
(ICR F) source J1317+34 (z = 1.050; iHewitt fc Burbided 
1989) were alternately observed through each eight-hour 
long VLBI run. The observations consisted of ~ 145 s scans 
on SN 2001gd and ~ 95 s scans on J1317+34, plus ~ 60 ad- 
ditional seconds of antenna slew time to make a duty cycle 
of ~300 s. In each observing run, the total on-source time 
was Ri4.2h and Ri2.5h for SN 2001gd and J1317+34, respec- 
tively. J1317+34 was also used as the amplitude calibrator 
for SN2001gd. The sources 4C 39.25 and J1310+32 (both 
ICRF sources) were observed as fringe finders. 

The correlated data were analyzed using the Astronom- 
ical Image Processing System (AIPS). The visibility ampli- 
tudes were calibrated using the system temperature and gain 



information provided for each telescope. The instrumental 
phase and delay offsets among the 8 MHz baseband convert- 
ers in each antenna were corrected using a phase calibration 
determined from observations of 4C 39.25 and J1310+32, for 
our VLBI observations on 26 June 2002 and 8 April 2003, re- 
spectively. We fringe-fit the data for the calibrator J1317+34 
in a standard manner; these data were then exported to the 
Caltech program DIFMAP (Shepherd et al. 1995) for imag- 
ing purposes. The final source image obtained for J1317+34 
was then included as an input image in a new round of fringe- 
fitting for J1317+34 within AIPS. In this way, the results 
obtained for the phase delays and delay-rates for J1317+34 
were nearly structure-free. These new values were then in- 
terpolated and applied to the source SN 2001gd using AIPS 
standard procedures. The SN 2001gd data were then trans- 
ferred into the DIFMAP program. We used standard phase 
self-calibration techniques using a time interval of 3 minutes, 
to obtain the images shown in Sect. [J] 



2.2 VLA measurements 

Our second observing epoch (2003.27) included the phased- 
VLA (in D-configuration) as an element of our VLBI array. 
We allocated about one hour of the VLA time to determine 
the radio spectrum of SN 2001gd. We observed at 1.4, 4.9, 
8.4, 15.0, 22.5, and 43.3 GHz, in standard continuum mode. 
We observed in both senses of circular polarization and each 
frequency band was split into two intermediate frequencies 
(IFs), of 50 MHz bandwidth each. We used 3C 286 as the 
primary flux calibrator (assumed of constant flux density) 
at all frequencies. We observed SN 2001gd phase-referenced 
to J1310+32 (a nearby VLA phase-calibrator source) at all 
bands, except at 8.4 GHz, for which we used our VLBI am- 
plitude calibrator, J1317+34, as the phase-reference source. 
J1310+32 and J1317+34 served also as the phase- reference 
sources for the system calibration. The flux densities re- 
ported here were obtained by combining the data from both 
IFs at each frequency band. As with our VLBI data, we 
used standard calibration and hybrid mapping techniques 
to obtain a VLA image of SN 2001gd and its host galaxy, 
NGC 5033, at the observing frequency of 8.4 GHz (see 
Fig. ID- 



2.3 Two-epoch XMM-Newton X-ray observations 

NGC 5033 was observed with the XMM-Newton Observa- 
tory on 2 July 2001 (PI: Turner, Obs. ID: 0112551301) and 
on 18 December 2002 (PI: Quataert, Obs. ID: 0094360501). 
Both observations used th e EPIC/MOS1, EPIC/MOS2, and 
EPIC /pn CCD cameras dStriider et al.ll200lt iTurner et al] 
2001); the two EPIC/MOS cameras were operated in the 
Prime Full- Window Mode, while the EPIC/pn camera was 
operated in the Extended Prime Full- Window Mode. The 
useful exposure times were ~7.6 ks and ~11.6 ks for the 
EPIC/MOS cameras and ~4.0 ks and ~10.0 ks for the 
EPIC/pn camera for the observations ID 0112551301 and 
0094360501, respectively. The observation ID 0112551301 
used the Thin Filter for the EPIC/MOS observations; the 
Medium Filter was used for all the other observations. 

We retrieved the XMM-Newton pipeline products from 



High-resolution observations of SN 2001gd in NGC 5033 3 



g 37 
< 



8.4 GHz 



— O SN 2001gd 




NGC 5033 



30 25 20 

RIGHT ASCENSION (J2000) 



Figure 1. Hybrid image of the galaxy NGC 5033 and its su- 
pernova SN 2001gd made at 8.4 GHz with the Very Large Ar- 
ray (VLA), from observations on 8 April 2003. The contours are 
(3,3 V3,9,...)x 16 fijy beam -1 , the off- source rms flux density per 
beam. The peak of brightness of the image corresponds to the nu- 
cleus of NGC 5033 and is 2.44 mjy beam -1 . The supernova is the 
bright point-like source northwards of the nucleus of NGC 5033. 
The principal dimensions (full width at half maximum) of the 
restoring beam are 10" 1 x8'.'7, with the beam's major axis ori- 
ented along a position angle of -61° (see inset in the lower left hand 
corner). 1" in the image corresponds to a linear size of 64 pc. 
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Figure 2. Contours of soft 0.25-2.5 keV X-rays from SN 2001gd 
overlaid on a Second-Epoch Digitized Sky Survey (DSS-2) Blue 
plate of NGC 5033. The X-ray image was extracted from the 
18 December 2002 XMM-Newton observation of NGC 5033 (Obs. 
ID 0094360501). The EPIC/MOS and EPIC/pn images have been 
combined to obtain a higher S/N image. The raw EPIC image has 
been further adaptively smoothed using Gaussian profiles with 
FWHM ranging from 1"5 to 9". The three contours in the image 
for SN 2001gd are drawn at (3, 8, and 13) X 0.15 cnts/arcsec 2 , the 
image off-source rms. 



the XMM-Newton Science Archive 1 . Initial processing and 
analysis were performed using the XMM-Newton Science 
Analysis Software (SAS ver. 6.0.0) and the calibration files 
from the Calibration Access Layer available on 13 April 
2004. Subsequent spectral analysis was carried out with 
XSPEC. 

We examined the light curve of the EPIC/MOS and 
EPIC/pn in the energy range above 10 keV to assess the ef- 
fect of any high energy particle flux on the background level. 
Although the light curve of the observation ID 0094360501 
on 18 December 2002 does not show any period of enhanced 
background, the observation ID 0112551301 on 2 July 2001 
is severely affected, with a background level 40 times higher 
than on the 18 December 2002 observation. Finally, the data 
were filtered to remove poor event grades (most likely spu- 
rious X-ray events) . 

We extracted images from the EPIC/MOS and 
EPIC/pn observations in different energy bands with pixel 
size of 1". The images extracted from the observation on 2 
July 2001 do not show any X-ray source at the location of 
SN 2001gd. On the other hand, the image in the 0.25-2.5 
keV band extracted from the observation on 18 December 
2002 (Fig. |5J does show a compact, soft X-ray source at 

1 The XMM-Newton Science Archive is supported by ESA and 
can be accessed at http://xmm.vilspa.esa.es 



the location of SN 2001gd, with 19±5 counts in the 0.3-2.0 
keV band. In the energy band above 2 keV, there are only 
2±2 counts, which corresponds to a hardness ratio of 0.11° '\ 
between the energy band above 2 keV and the total XMM- 
Newton energy band. The facts that the pre-supernova ob- 
servation yielded a non-detection and the post-supernova 
observation yielded a detection, coupled with the spatial 
coincidence and soft X-ray spectrum imply that the X-ray 
source is SN 2001gd, and that the supernova explosion likely 
happened after 2 July 2001. 



3 VLBI IMAGES OF SN 2001GD 

We used phase-reference techniques at both VLBI epochs to 
obtain images of SN 2001gd. The use of phase-referencing 
allows detection and imaging of very faint sources like su- 
pernovae, as it effecti vely increases the coheren ce time from 
minutes to hours (e.g.. lBeaslev fc Conwavll995l) . In our case, 
we made phase-reference observations of our target source, 
SN 2001gd, with respect to the relatively strong, nearby 
QSO J1317+34, which is 2°22' distant. 

We display in Fig. our 8.4 GHz VLBI image of 
J1317+34 on 26 June 2002. The source shows a one-sided 
core-jet structure at milliarcsecond (mas) scales, with the 
jet extending northwards to ~ 9 mas. As mentioned in 
Sect. 12.11 we subtracted at each epoch the phase contri- 
bution of J1317+34 (due to its extended structure) from 
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Figure 3. 8.4 GHz Very-Long-Baseline Interferometry (VLBI) 
image of J1317+34 on 26 June 2002. The contours are 
(3,3%/3,9,...)x80/iJy beam -1 , the off- source rms flux density per 
beam. The peak of brightness of the image is 184 mjy beam -1 
and the restoring beam (bottom left in the image) is 1.29 mas 
xO.61 mas at a position angle of -14?3. The total flux density 
recovered in the image is ~ 0.33 Jy. North is up and east is left. 
One milliarcsecond in the image corresponds to a linear size of 
8.5 pc. 



the fringe phase, and used the position of its phase-center 
(essentially coincident with the peak of the brightness distri- 
bution), as our reference point for phase-reference imaging. 

Figure 2] shows our VLBI images of SN 2001gd, which 
showed no structure at both epochs. The coordinates of 
the supernova are offs et from those provided by the VLA 
iStockdale et alJ 120021) . and used as a priori VLBI coor- 
dinates in the correlator, for the supernova position, by a 
mere ~3 mas in right ascension and ~1 mas in declination. 
Since the standard er ror of the VLA position provided by 
IStockdale et alJ ll2002l) was 0.2 arcsec in each coordinate, the 
(small) differences between the VLBI and VLA determined 
coordinates would be a highly improbable event. It is more 
likely that the VLA uncertainties quoted bv lStockdale et al] 
(2002) were far too conservat ive. Indeed, a reanalysi s of the 
8.4 GHz VLA observations of lStockdale et all i2002i) results 



in 1 a uncertainties of 0.015-0.0020 arcsec (C. Stockdale, 
private communication) for the VLA position of SN 2001gd, 
a bout a factor of ten sm aller than the uncertainties quoted 
in lStockdale et alJ d2002l) . 



4 RESULTS AND DISCUSSION 

4.1 Radio light curve and spectral behaviour 

Simultaneous multiwavelength radio observations of super- 
novae are important to help characterize the relevant phys- 
ical processes taking place. Since we were granted VLA ob- 
serving time for our second VLBI epoch on 8 April 2003 
we were able to characterize the radio spectral energy dis- 
tribution of SN 2001gd. The flux density measurement er- 
ror given for the VLA measurements in Table represents 
one statistical standard deviation, and is a combination 
of the off-source rms in the image and a fractional error, 
e, included to account for the normal inaccuracy of VLA 
flux density calibration and possible deviations of the pri- 
mary calibrator from an absolute flux density scale: the 
final errors, cry, as lis ted in Table ^ are ta ken as a* = 
(eS ) 2 + oo (see, e.g., IStockdale et al.ll2003l) . where So is 
the measured flux density, ao is the off-source rms at a 
given frequency, and e=0.10 at 1.4 GHz, 0.05 at 4.9 GHz, 
0.05 at 8.4 GHz, 0.075 at 15.0 GHz, and 0.10 at 22.5 and 
43.3 GHz. We fit our data using an implementation of the 
nonlinear weighted-least-squares Marquardt-Levenberg al- 
gorithm. We plot in Fig. two model fits to the data in 
Tabled The solid line corresponds to a fit to the data with 
a power-law spectrum that is partially suppressed by free- 
free absorption from a homogeneous screen of ionised gas 
[S v = Si v a exp( — Tff,„)], where Si is the flux density at 
1 GHz, in mjy; v is the observing frequency, in GHz; and 
TS,w = Av~ 2A is the free-free optical depth at the ob- 
serving frequency v. The dashed-dotted line corresponds 
to a pure synchrotron spectrum fit. For our first model, 
we obtain the following weighted- least-squares ("best-fit") 
values for the parameters: £i=(9.1±1.6) mjy, a=-1.0±0.1, 
A=(1.0±0.3) GHz 2 ' 1 . (The quoted uncertainties correspond 
to lcr.) The observed spectral index is typical of the opti- 
cally thin phase of Type II supernovae. The corresponding 
optical depth at 8.4 GHz is r ffjS . 4 =(l.l±0.3) x 10" 2 . This es- 
timate can be used to constrain the electron temperature, 
T e =10 T5 K, and mass-loss rate parameter of the supernova 
progenitor, M/v w . (Each letter symbol with a subscript in- 
dicates the value of the corresponding quantity in units of 10 
to the power given by the subscript; for example, T5 above 
denotes T e in units of 10 5 K.) Indeed, for a steady presu- 
pernova mass loss rate, M=10 _5 M_5 MQyr -1 , and wind 
speed, i) m =10i)io km s -1 , the free- free absorption optical 
depth of the unshocked ionised gas is llLundavist fc Franssonl 

Ts,u « 0.17 g s fiT- 3/2 v~ 2 M! 5 v^ 2 r^ 3 (1) 

where v is the frequency in GHz, r — 10 16 rig cm, is the 
radius of the circumstellar shock, gg is the free-free Gaunt 
factor, and fj, = [1 + 2n(He)/n(H)]/[l + 4n(He)/n(H)] is 
the mean molecular weight for the presupernova wind. For 
solar abundances, n(He)/n(H)=0.1 and fi— 0.86, while for 
n(He)/n(H) = l, as may have been the case for SN 1993J 
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Figure 4. Top: Very-Long-Baseline Interferometry (VLBI) im- 
age of SN 2001gd on 26 June 2002. Contours are drawn at 
(3,3 v / 3,9,...)xlOAtJy beam -1 , the off- source rms. The restoring 
beam (top right) is 1.41 X 0.48 mas 2 at a position angle of -13? 4. 
The peak of brightness is 3.64 mjy beam -1 , and the total flux 
density recovered for the supernova is 3.83 mjy. Bottom: Very- 
Long-Bascline Interferometry (VLBI) image of SN 2001gd on 8 
April 2003. Contours are drawn at (3,3 V3,9,...) X 11 [ijy beam -1 , 
the off-source rms. The restoring beam (top right) is 1.55 X 0.59 
mas 2 at a position angle of -23? 7. The peak of brightness is 
944 /ijy beam -1 , and the flux density recovered for the super- 
nova is 1.02 mjy. In both panels, north is up and east is left; the 
origins are taken from the VLA observations. 
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Figure 5. Model fits to the radio spectrum of SN 2001gd on 8 
April 2003. The solid line corresponds to a synchrotron spectrum 
fit, partially suppressed by free-free absorption, whilst the dashed- 
dotted line corresponds to a pure synchrotron spectrum fit for all 
frequencies, except 1.4 GHz. The error bars denote ±ler values. 



jBaron et alJll994l) . /x=0.6. Using n 6 = 3.6 (see Sect. IO 
below) and evaluating gg at 8.4 GHz, we obtain T$ tV 



2.1 X 10 -4 (m/0.6) T~ 



2 . iFransson et alJ (Il99' 



and IFransson fc Biornssonl 1^)98) showed that Tg sa 2 and 
(M-s/vio) ~5 were needed to explain the X-ray and ra- 
dio emission of SN 1993J. Using the above values of 
and (M- 6 /vi ) for SN 2001gd on 8 April 2003, we obtain 
from Equation Q Tfi,8.4 ~ 2 x 10~ 3 , which is a factor of 
three too low to fit the lower end of our value for ts,&a- 
We explored the space of possibilities within 2 a of our best- 
fit values for the parameters Si, a, and A (see above). We 
found very good matches for ts : sa for the pairs of values 
(T e =3 X 10 4 K; M=2.5 x 10 -5 M yr -1 ) and (T e =2 X 10 5 K; 
M=l x 10" 4 M©yr -1 ). Unfortunately, the coupling of T e 
and (M/v w ) prevents us from favouring one pair of val- 
ues over any other solely based on the radio observations, 
so we limit ourselves to pointing out that the values of 
T e and M for SN 2001gd are likely to be in the range 
T e =(3 - 20) x 10 4 K, M=(2 - 10) x 10" 5 M Q yr -1 . As we 
will see in Sect. 14.21 the available X-ray data favour models 
with M < 5 x 10 -5 MQyr -1 , which would imply electron 
temperatures of T e < 7 x 10 4 K to satisfy our radio obser- 
vations. 



4.2 X-ray emission from SN 2001gd 

Mode ls for the circumstella r interaction around supernovae 
(e.g., IFransson et al.l Il996l) predict the existence of two 
components in X-rays: a soft component associated with 
the reverse shock of the supernova, T IS , and a hard com- 
ponent ass ociated with the circum stellar shock, T cs , e.g., 
SN 1995N jMucciarelli et alJl2004|) . The 18 December 2002 
(i-to=472 d) XMM-Newton observations of NGC 5033 can 
be used to investigate the circumstellar interaction around 
SN 2001gd and to estimate its X-ray luminosity. However, 
the spectrum has too small a number of counts to allow 
spectral fitting. Instead, we compared the count rate and 
hardness ratio of the observed XMM-Newton EPIC / pn spec- 
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Table 1. Results from VLA observations on 8 April 2003 



Source Name 






Flux density (m Jy) 








43.3 GHz 


22.5 GHz 


15.0 GHz 8.4 GHz 


4.9 GHz 


1.4 GHz 


SN 2001gd 
1310 + 323M 
J1317 + 34( 2 > 


£ 0.30+ 
2584 ± 258 


0.39 ±0.12 
2911 ±291 


0.63 ±0.13 1.02 ±0.05 
2397 ± 178 

384 ± 19 


1.85 ±0.24 
1656 ± 83 


3.89 ±0.40 
1535 ± 154 



W Phase and secondary flux density calibrator at all frequencies, except at 8.4 GHz. ( 2 ) Phase and secondary flux density calibrator 
at 8.4 GHz. t The quoted flux density at 43.3 GHz is an upper limit, corresponding to three times the off-source rms noise in the image. 
The quoted uncertainty for each flux density value corresponds to lc (see Sect. |4~T1 for details). 



20 



10 




20 



3 10 



o 

20 



10 



-i | i i i i | 
- (a) 



\ — III 



— i | i i 
s = 2 
n=10 
M_ 5 = 2.5 



mitipit, 



(b) 



~1 r I I 



s = 2 
n=10 
M_ 8 = 15 



(c) 





0.5 1 5 

Energy (keV) 



Figure 6. XMM-Newton EPIC / pn spectrum of SN 2001gd on 18 
December 2002, overplotted by thin plasm a X-ray emission mod - 
els for a supernova radiative reverse shock iFransson et al.ll996t) . 
The different sets of parameters are displayed in the correspond- 
ing panel. The uncertainties shown are one statistical standard 
deviation. The energy bin-width is 0.5 keV; therefore, 1 cnt in 
the spectrum corresponds to ~ 2.0 X 10 -4 cnt s _1 keV - . 



trum with these i mplied from the a nalyti cal models for X- 
ray supernovae bv iFransson et ail (I1996T) . Figure HJ shows 
the XMM-Newton spectrum of SN 2001gd on 18 Decem- 
ber 2002, overplotted by several optically thin plasma X-ray 
emission models. According to these models, both column 
density absorption and electron temperature are sensitive 
to the circumstellar density profile (p CS m ~ r~ s ), ejecta den- 
sity profile (p c j ~ r~ n ), shock speed (u s h = 10 4 va km s _1 ), 
and mass-loss rate (Af=10~ 5 M_s M yr -1 ). As for the ex- 
pansion velocity, H a spectra taken on December 2001 and 
January 2002, indicate whwhm » 7800 ± 1000 km s _1 (T. 
Matheson, private communication). Considering a likely su- 
pernova deceleration from this date to the time when the 
X-ray observations were obtained, we adopted a supernova 
shock speed of V4=0.6 on 18 December 2002. We also kept 
fixed the presupernova wind speed at Wio=l, while we varied 
s, n, and M to investigate their effects on the X-ray emis- 
sion. W e note that the analytical models of IFransson et all 
( 1996) do not predict the exact amount of X-ray luminosity, 
but the reverse shock total luminosity, L rs . A close inspec- 
tion of their results indicates that the X-ray luminosity is 
~20% of L ra , and thus we adopted an X-ray luminosity of 
this fraction of L rs to estimate the XMM-Newton EPIC/pn 
count rate. 

The supernova model that shows hardness ratio and X- 
ray luminosity most consistent with the XMM-Newton spec- 
trum of SN 2001gd on 18 December 2002 has kT IS » 1.1 
keV, s=2, n=10, and Af_s=2.5. For illustration purposes, 
we compare in Fig. the spectral shape of this model (top 
panel) with the observed X-ray spectrum. The above pa- 
rameters agree well wit h what is expected fo r a supernova 
radiative reverse shock (IFransson et al]ll996l) . Our value of 
M_5=2.5 for the mass-loss rate of SN 2001gd from the 
X-ray observations may be compar e d wit h the value of 
M_5=3 obtained by IStockdale et al.l J2003) from their ra- 
dio light curve modelling, and with our own estimates in 
Sect. 14.11 The corresponding column density of the cool- 
ing shell between the reverse shock and the contact surface 
is N coo \ = 5.5 x 10 21 cm~ 2 , and that of the circumstellar 
shock is N cs = 1.6 x 10 21 cm -2 . We note that the supernova 
shock speed, v B h, and the ejecta density profile, n, are cou- 
pled, so that different combinations within a certain range 
of these parameters reproduce the observed hardness ra- 
tio and count rate. Models with supernova shock speeds as 
high as Vi ~ 1.2, would require n « 14 to reproduce the 
observed hardness ratio, but result in unreasonably high 
XMM-Newton EPIC/pn count rates. On the other hand, 
models with supernova shock speeds as low as v± m 0.2, 
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would require n > 5 to reproduce the observed hardness 
ratio, but result in too low XMM-Newton EPIC/pn count 
rates. We also explored models with mass loss rates and cir- 
cumstellar density profiles in the ranges 0.1 M-5 ^ 30 
and 1.1 < s < 2.9, respectively. Models with M_ 5 < 0.5, 
while predicting consistent hardness ratios, are unreason- 
able as their cooling times are much greater than the age of 
SN 2001gd. Models with s > 2.4 predict significantly lower 
XMM-Newton EPIC/pn count rates than observed. Models 
with M_5 > 15 (middle panel in Fig.[SJ or s < 1.4 (bottom 
panel in Fig. |SJ have hardness ratios between the energy 
band above 2 keV and the total XMM-Newton energy band 
that are > 0.6, well above the observed hardness ratio of 
O.llo.i (see Sect. ED. Therefore, models with M_ 5 < 0.5 
or M_5 > 15, and models with s < 1.4 or s > 2.4 seem less 
likely to describe the X-ray emission from SN 2001gd than 
the model shown in the top panel of Fig. [Ij] 

At the adopted distance of 13.1 Mpc, the unabsorbed 
flux from SN 2001gd is F x = (7 ± 2) x 10 -14 ergs cm -2 s _1 
in the 0.3-5 keV band, which corresponds to a luminos- 
ity of L x = 47rD 2 fx= (1.4±0.4) x 10 39 L> 2 3 ergs s _1 , where 
L>=13.1 D13 Mpc is the distance to the supernova. The above 
results are in agreement with expecta tions from X-ray emis - 
sion from Type II supernovae (e.g., iFransson et al.lll996t) . 
and seem to favour a relatively shallow ejecta density profile 
and a mass-loss rate about a factor of two lower than for 
SN 1993J. 

The 2 July 2001 XMM-Newton observations of 
SN 2001gd did not yield detectable X-ray emission, but can 
be used to derive an upper limit of the isotropic X-ray lu- 
minosity of SN 2001gd at that time. The 3cr upper limit on 
this luminosity is Lx = 7.4 x 10 38 erg s _1 . Note, however, 
that this upper limit to the X-ray luminosity was obtained 
with a significantly shorter exposure time, and with a much 
higher background X-ray emission level, compared to that 
for the observations on 18 December 2002. 

Finally, the expected monochromatic flux density 
at 4.9 GHz on 18 December 2002 [t - t =472 d) was 
w (2.5±0.3) mJy. Using a = -1.0 from our best-fit VLA 
spectrum on 8 April 2003, we obtain an isotropic radio lu- 
minosity of Lr = (8 ± 2) x 10 36 D 2 3 ergs s _1 between 1.4 
and 43 GHz. The values of L x and L R put SN 2001gd close 
to SN 1993J iChevalierll2003l) regarding radio and X-ray lu- 
minosities, in agreement with the fact that SN 2001gd dis- 
played a very similar behaviour to SN 1993J in the optical 
band. 



4.3 Expansion of SN 2001gd 

The VLBI images of SN 2001gd at 8.4 GHz (Fig. H show 
a compact, unresolved radio emitting structure. Therefore, 
we resorted to model fitting to estimate the supernova an- 
gular diameter, 9, at each epoch. For this purpose, we 
used the model fitting procedure included in the Caltech 
DIFMAP package, which fits a given model directly to the 
real and imaginary parts of the observed interferometric 
visibilities, using the Levenberg-Marquardt non-linear least 
squares minimi zation technique . Inspir ed by the findings for 
SN 1993J (e.g.. iMarcaide et alJll995albL Il997t iBartel et alJ 
2000), we assumed spherical symmetry in our model fits, 
and considered the following models: (i) an optically thin 



Table 2. Angular diameter estimates of SN 2001gd 



Model 


Angular diameter (/t as) 




2002.48 


2003.27 


OPTICALLY THIN SPHERE 


390 ± 80 


440 ± 100 


UNIFORMLY BRIGHT DISK 


350 ± 60 


410 ± 80 


OPTICALLY THIN SHELL 1 " 


330 ± 60 


370 ± 80 



t The shell width is 25% the outer radius. 

sphere; (ii) a uniformly bright, circular disk; and (iii) an op- 
tically thin she ll of width 25 % its outer shell radius (e.g. 
iMarscheJllOSfih . Table 2 summarizes our results, where the 
uncertainties (la) represent one statistical standard devia- 
tion plus an estimate of modelling errors. All three models 
give similarly good fits to the (u, v) data (x? cd = 0.8-0.9), 
so in principle we cannot rule out any of them. Our angular 
diameter estimates nominally suggest a relatively strong de- 
celeration in the expansion of SN 2001gd, but values of the 
deceleration parameter m (8 tx t m ) between and 1 are all 
within the uncertainties, and therefore we cannot draw any 
conclusion from our data. 

4.4 Energy budget and magnetic field in 
SN 2001gd 

The high brightness temperatures inferred from our VLBI 
observations for SN 2001gd [(5.3±1.2)x 10 8 K on 26 June 
2002, and (1.1T0.3) x 10 s K on 8 April 2003], indicate a 
non-thermal, synchrotron origin for the radio emission from 
SN 2001gd. We can estimate a minimum total energy in rel- 
ativistic particles and fields, and an equipartition magnetic 
field for SN 2001gd, by assuming equipartition between fields 
and particles. The mini mum total energy and its equiparti- 
tion magnetic field are dPacholczvklll970l) : 

£ m in = C 13 (l+^ 4/7 3/7 i? 9/7 4 /7 (2) 

B min = (4.5 c 12 /0) 2/r (1 + V) 2/? R~ 6/7 4 /? (3) 

where Lr is the radio luminosity of the source; R, its linear 
radius; C1 2 and C13, slowly- varying functions of the spectral 
index, a iPacholczvkll970T) : cj>, the filling factor of fields and 
particles; and ip the ratio of the (total) heavy particle energy 
to the electron energy. This ratio depends on the mechanism 
that generates the relativistic electrons and ranges from ip w 
1 to ip = m p /m e » 2000, where m p and m e are the proton 
and electron mass, respectively. 

From our VLA observations on 8 April 2003, we de- 
termined a spectral index a = — 1.0 ± 0.1 (S„ oc v a \ 
see Sect. 14. in between 1.4 and 43 GHz, and Ss.4GHz = 
1.02 ± 0.05 mJy. With these values, we obtain Lr = (6.0 ± 
0.3) x 10 36 D 13 ergs s" 1 between 1.4 and 43 GHz. We used 
0=0.6 0o.6 in our calculations, which corresponds approxi- 
mately to the filling factor for the shell model in Table 2. In 
any case, the estimates of L m i n and L? m i n are relatively insen- 
sitive to the precise value of <f>. The radius of SN 2001gd on 8 
April 2003 is, for the shell model, R = 3.6 x 10 16 L>i 3 8i$ 5 cm, 
where 8iss is the angular radius of the supernova in units of 
185 mas. With these values, we get from Eq. © and (|3J 

L min » 1.8 x 10 46 (1 + ^) 4 ^ 7 3/ 6 7 8% 7 5 DZ /7 ergs 
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B min « 40 (1 + ^) 2/r <f«T QisV D~ 3 2/7 mG 

Since 1 < tj) < 2000, E mln is in the range (2.6 x 10 46 - 1.4 x 
10 48 )erg, and the equipartition magnetic field in the range 
(50-350) mG. Thus, for an energy of 10 51 ergs for the ex- 
plosion of SN 2001gd, the fraction of energy necessary to 
power the radio emission of SN 2001gd is quite modest. 
The upper range of the equipartition magnetic field is lower 
than, but not far from, the lowest part of the range of val- 
ues obtained for SN 1993J at similar radii (B ~ 1.7 G, 
Fra nsson fc Biornssonlll998l : B G. IPerez- Torres et alJ 

Since the kinetic energy density in the wind is 
likely larger than the circumstellar magnetic field en- 
ergy density (e.g., iFransson fc Biornssonl Il998f) . we have 
pcsm«w/2 <; B 2 Brn /8n, where p csm and B csm are typical of 
the density and magnetic field in the circumstellar medium, 
respectively. By assuming a standard wind density profile 
(p oc r -2 ), we can write the above expression as 

B csm < (Mv w ) 1/2 r- 1 » 2.5 (M_ 5 v 10 ) 1/2 r^ 1 mG (4) 

where r\§ = r/10 16 cm. For rw = 3.6 and M_s = 2.5, 
we obtain B csln < 1 mG, which is a factor about 50 to 350 
times smaller than the equipartition field, and shows that 
the magnetic field inferred for SN 2001gd cannot originate 
solely by compression of the existing circumstellar magnetic 
field, which would increase the field only by a factor of four 
(e.g., iDvson fc Williams! Il980l) Large amplification factors 
of the magnetic field have also been found for other ra- 
dio S Ne, e.g., SN1993J (amplification factors of a few hun - 
dred ; IFransson fc Biornssonll 19981 . IPerez- Torres et alj|200lfl , 
or SN 1986J (amp lification factors in the range 40—300; 
IPerez- Torres et alJ [20021. Thus, if equipartition between 
fields and particles exists, amplification mechanisms other 
than compression of the circumstellar ma gnetic field need 
to be invoked, e.g., turbule nt amplification JChevalierlliiil 
IChevalier fc BlondinllT995T) . to explain the radio emission 
from supernovae. 



5 SUMMARY 

We have presented the results of the first two 8.4 GHz very- 
long-baseline interferometry observations of SN 2001gd in 
NGC 5033, complemented with VLA and XMM-Newton ob- 
servations. We summarize our main results as follows: 

• The radio structure of SN 2001gd is not resolved, either 
on 2002.48 (t - to = 295 days; i =3 September 2001) or on 
2003.27 (t — to = 582 days). We used the interferometric vis- 
ibility data to estimate the angular sizes for the supernova, 
as well as constraints on its expansion speed from optical 
measurements. While our data nominally suggest the pos- 
sibility of a relatively strong deceleration in the expansion 
of SN 2001gd (m < 1, oc i m ), solutions with m = 1 are 
also permitted, and therefore we cannot draw any conclusion 
from these data. 

• Our VLA data on 8 April 2003 can be well fit by an opti- 
cally thin, synchrotron spectrum (a = — 1.0±0.1; S v oc u a ), 
which is partially absorbed by thermal electrons. The radio 
spectrum and the light curve at 1.4 GHz indicate that the 
synchrotron turnover frequency is near 1 GHz. We obtain a 
supernova flux density of (1.02L0.05) mjy at the observing 



frequency of 8.4 GHz. At an adopted distance of 13.1 Mpc to 
NGC 5033, this flux density implies an isotropic radio lumi- 
nosity of (6.0±0.3)x 10 36 ergs s" 1 between 1.4 and 43 GHz. 
We also used our best-fit to the VLA radio spectrum to infer 
the most likely ranges for the electron temperature in the 
wind, T e , and the mass-loss rate, M. From our VLA obser- 
vations on 8 April 2003, we find T e = (3 - 20) x 10 4 K, and 
M=(2 - 10) x 10" 5 M©yr _1 . (Note the constraint on the 
combination of these values; see Sect. l4.l1 and l4.2T ) 

• By assuming equipartition between fields and particles, 
we estimate a minimum total energy in relativistic particles 
and magnetic fields in the supernova shell of -E m i n «(0.3- 
14) x 10 47 erg, which corresponds to an equipartion aver- 
age magnetic field there of B m i n ~ (50-350) mG. We find 
that the average magnetic field in the circumstellar wind of 
SN 2001gd is Bcsm ^ 1 mG at a radius from the center of 
the supernova explosion of » 3.6 x 10 16 cm. Since compres- 
sion of this existing magnetic field by the supernova shock 
front could enhance it by only a factor of four, powerful am- 
plification mechanisms must then be acting in SN 2001gd, 
to account for the magnetic fields responsible for the syn- 
chrotron radio emission. 

• We used XMM-Newton archival data to estimate the 
X-ray flux of SN 2001gd on 2 July 2001 and 18 December 
2002. The data from the X-ray observations on 2 July 2001 
are below the noise level, while we detected SN 2001gd on 18 
December 2002, with a flux of (7 ± 2) x 10" 14 erg cm" 2 s" 1 
in the 0.3-5 keV band. The emission is chiefly from a soft 
component (<2 keV), associated with the reverse shock of 
the supernova. The corresponding isotropic X-ray luminos- 
ity is L x = (1.4±0.4)xl0 39 erg s" 1 in the 0.3-5 keV band. 
The X-ray spectrum is consistent with expec tations from X- 
ray e mission from Type II supernovae (e.g.. IFransson et alJ 
1996), and suggests the interaction of a relatively shal- 
low supernova ejecta density profile (p e j ~ r~ n ;n ~ 10), 
with a standard circumstellar wind density profile (p csm ~ 
r~ s ; s ~ 2.0), characterized by a presupernova mass-loss rate 
of M « 2.5 x 10~ 5 M yr" 1 for SN 2001gd. 

SN 2001gd resembles SN 1993J in its radio and X-ray 
emission, indicating that a similar circumstellar interaction 
is taking place. Since available radio supernovae for VLBI 
studies are so rare, it is important that the few of them 
that allow such studies be monitored if we are to better un- 
derstand the radio supernova phenomenon. Further VLBI 
observing epochs of SN 2001gd are necessary to trace the ex- 
pansion history of the supernova, and to determine whether 
SN 2001gd is decelerating. However, the task is most chal- 
lenging, as SN 2001gd appears to fade very quickly in radio. 
We estimate that SN 2001gd will have a flux density at 8.4 
GHz of about 0.3-0.4 mjy around March- April 2005. There- 
fore, the chances of a successful radio imaging of SN 2001gd 
are small, and will necessarily require the use of an array 
including the world's most sensitive antennas. 
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